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We propose a finite temperature Landau theory that describes competing orders and interlayer 
tunneling in cuprate superconductors as an important extension to a corresponding theory at zero 
temperature [Nature 428, 53 (2004)], where the superconducting transition temperature T c is defined 
in three possible ways as a function of the zero temperature order parameter. For given parameters, 
our theory determines T c without any ambiguity. In mono- and double-layer systems we discuss the 
relation between zero temperature order parameter and the associated transition temperature in 
the presence of competing orders, and draw a connection to the puzzling experimental fact that the 
pseudo-gap temperature is much higher than the corresponding energy scale near optimum doping. 
Applying the theory to multi-layer systems, we calculate the layer-number dependence of T c . In a 
reasonable parameter space the result turns out to be in agreement with experiments. 



The essential phenomenology of high-T c superconduc- 
tors are the dome-shaped superconducting phase diagram 
as a function of doping x, and the existence of a pseudo- 
gap normal state. The mechanism of the novel super- 
conductivity remains to be a challenge to researchers. 
The fact that superconductivity follows from doping anti- 
ferromagnetic Mott insulators and that the transition 
temperature depends on doping non-monotonically sug- 
gest the importance of the effects of strong correlation 
among the electrons, as emphasized firstly by Ander- 
son. [l| A theory in this line must go beyond the usual 
Bardeen-Cooper-Schrieffer mean field theory. However, 
there is also a possibility that the phenomenon is not 
that exotic, and a modification to the conventional the- 
ory by including a competing order may do the job. This 
is the attitude taken by Chakravarty, Kee and Volker 
0, who proposed a zero temperature Landau theory 
to explain the general phase diagram, and in particu- 
lar the copper-oxide layer-number N (within a unit cell) 
dependence of superconducting transition temperature 
T c of homologous series of cuprate superconductors. 
The new ingredients are just a competing order, e.g., 
the d-density-wave (DDW),Q and inter-layer tunnelling 
that enhances superconductivity. It is therefore im- 
portant to judge how robust the conclusions are against 
the weak points in the theory. Indeed, a few of them 
are debatable even within the mean field framework it- 
self. 1) When two orders compete with each other, it 
is not clear whether there is still a definite connection 
between the zero temperature order parameter and the 
transition temperature even if this is the case in the ab- 
sence of competition between the orders. 2) In the case 
of multi-layer systems, it is not clear which representa- 
tive of the modulated order parameters (due to charge 
redistribution) is most appropriate to be related to the 
transition temperature, even if one assumes that the zero 
temperature order parameter scales with the transition 
temperature. 3) Moreover, the underlying motivation for 
a competing order to superconductivity is the fact that 



the normal state pseudo-gap seems to be independent of 
the pairing gap. From the elaborate collection of data in 
Ref.[3 the pseudo-gap energy scale extrapolates to zero 
at a doping level of x = 0.19, while the pseudo-gap phe- 
nomena certainly exists at T > T c even at the same dop- 
ing. This already indicates that the zero temperature 
value of the competing order, assumed to be responsi- 
ble for the pseudo-gap, does not scale with the transition 
temperature for the pseudo-gap itself (if the pseudo-gap 
temperature crossover is a phase transition at all). As a 
compromise, the authors in Ref.Q defined a gap that is 
the root-sum-square of both order parameters, and used 
it to represent the pseudo-gap temperature. In general, 
there is no microscopic basis for this gap-combination. In 
a microscopic tight-binding model, the two gaps do not 
combine this way once the DDW bands is doped away 
from half filling. 4) Finally the weak points beyond the 
mean field theory is of course the ignorance of quantum 
or thermal fluctuation of the order parameters. 

The last point can only be addressed by going be- 
yond mean field theory, as in a recent spin wave analysis 
of the anisotropic XY- model, which reproduces nicely 
the experimental T C (N) in homologous series of HBCO 
superconductors. 3J In this paper, we address points l)-3) 
listed above in the mean field framework under the mo- 
tivation that the importance of a successful mean field 
theory should not be underestimated. In order to do 
so, we propose a finite temperature Landau theory, and 
determines the transition temperatures for both orders 
unambiguously. The structure of the rest of the paper is 
as follows. A finite temperature theory for the supercon- 
ducting order or the DDW order alone is first proposed. 
We then introduce the coupling between these two or- 
ders, and discuss how competing orders in the case of 
mono-layer and double-layer systems, where no charge 
redistribution occurs, modify the relations between zero 
temperature order parameter and the transition temper- 
ature. By comparing the theoretical phase diagram with 
experiments we fix the parameter introduced in the finite 



2 



temperature theory. We also comment on point 3) raised 
above. Finally we extend the theory to the multi-layer 
systems, and calculate the T C (N). We find that in a rea- 
sonable parameter space, the theoretical results are in 
agreement with experiment, Q and lend a support to the 
conjecture from the zero temperature Landau theory. 
Let us begin with the simple case that the supercon- 
ducting order is the only order parameter. It is known 
that 2tp /T c = 3.52 in the weak coupling s-wave super- 
conductor, where ipQ is the zero temperature energy gap. 
For the case of weak coupling d-wave superconductor, 
this ratio is roughly 4.16. A Landau- type theory could be 
derived near the transition temperature from an effective 
microscopic model, but is a daunting task to access the 
zero temperature limit. We would therefore take a phe- 
nomenological attitude. We demand that the theory is 
of the form near T c , and while extended to lower temper- 
atures should give a qualitative temperature dependence 
of the order parameter, and in particular, should yield 
a prescribed ratio of ipo/T c . Without loss of generality, 
we re-scale the temperature so that ipo/T c = 1. Under 
these conditions, a suitable finite temperature Landau 
free energy density is as follows, 

f s = [a(x)+/3T 2 M 2 + ^\\ (1) 

Henceforth we use arbitrary units, and borrow the pa- 
rameterization in Ref. 2], with a(x) — 10(x— 0.3) and f3 = 
2. This simple model would yield T c = ipo = y—a(x) / (5. 
We note that near T c the coefficient of |^| 2 -term can be 
approximated as 2(3T C (T — T c ), which is of the desired 
form in a usual Landau theory. 

The same consideration can be applied for the DDW 
order alone, with the free energy density 

f D = [a{x)+qbT 2 ]D 2 + h -D i , (2) 

where a(x) — 27 '(x — 0.22), 6 = 2, and q is a new phe- 
nomenological constant. The transition temperature for 
the DDW order alone is Td = yj—a(x)/qb = D /^/q, 
where Dq is the bare zero temperature DDW order. Since 
we have taken the latitude to set tpo/T c = 1, there is no 
further freedom to set q = 1. In fact the value of q is 
given by q = (D /To) 2 / '(^>o /T c ) 2 independently of the 
temperature re-scaling. It is therefore understood that 
a smaller value of q means that the DDW order is more 
tolerant to thermal suppression, and vice versa. We take 
q as a new phenomenological parameter, which will be 
fixed by arguments below. 

Let us now couple the two order parameters to form a 
theory describing mono-layer systems, 

/= [ a (x) + PT 2 M 2 + ^ 

+ [a(x) + qbT 2 }D 2 + ^D i + gD 2 \4>\ 2 , (3) 



where g — 1.2 is the coupling constant. As has been il- 
lustrated in Ref.0, this theory yields zero temperature 
order parameters ipo(x) and D (x) that depends on the 
doping level x in very much the same way as the em- 
pirical T c (x) and T]j(x), provided that one assumes that 
Td oc y/tpQ + Dq. Here Td represents the pseudo-gap 
temperature scale. However, as has been pointed out, 
there is no microscopic basis for the gap-combination, 
and the relation between such a combined gap and the 
pseudo-gap temperature is not clear a priori. In our case, 
the transition temperature T c for i\) is determined by 

a(x) + (3T 2 + gD 2 = 0, (4) 

where D = y/-[a(x) + qbT^jb is the value of D at T = 
T c . Similarly, the transition temperature To for D is 
given by 

a(x) + bqT% + gi? = 0, (5) 

where ip = y/—[a(x) + f3T D ]//3 is the value of tp at 
T = To- Explicit expressions for T c and Td are not 
provided because D > and tfj > have to be im- 
posed in the above equations. We present T c (thick solid 
lines) and Td (thin solid lines) for typical cases, with 
q = 0.1,0.6, 1, 1.5 in Figs.l(a)-(d) respectively. We have 
also shown ipo (thick dotted lines) and Dq (thin dotted 
lines) in each figure for comparison to the transition tem- 
peratures. We observe that ipo/T c = 1 for all doping lev- 
els in Fig. 1(b). However, this does not hold in the other 
cases. Moreover, in the same case D /Td is apparently 
not a constant. The reason why Fig. 1(b) is a special case 
is because the parameters satisfy qb — g = 0, which gives 
rise to a line of critical temperatures for DDW at the 
zero-temperature critical doping level. For q < 0.6, as 
in Fig. 1(a), we observe a lobe-shaped DDW phase tran- 
sition line near x = 0.198, namely, two Tq's at a given 
doping. This is best understood in the limit of q = 0. 
Under this assumption DDW order is not suppressed by 
thermal effects but by the superconducting order. Near 
T c the latter is small and DDW may survive. However, 
with lowering temperature the superconducting order in- 
creases and eventually squeeze the DDW order. With a 
finite but small q the qualitative behavior is the same. 
Because such a feature is not observed yet, it is very un- 
likely that q < 0.6. For larger values of q as in Figs. 1(c) 
and (d), the resulting doping dependence of T c is very 
oblique, a feature that does not seem to appear in high-T c 
phenomenology. Basing on these judgements we believe 
a reasonable regime is given by 0.6 < q < 1. We empha- 
size an interesting feature in Fig. 1(b), where To survives 
above T c in the under-doped regime and drops abruptly 
to zero at x — 0.198, whereas Do vanishes nearby. The 
reason for this to occur is because D vanishes here not 
because of thermal suppression but the competing order 
tjj. This seems to account for the puzzling fact that even 
though pseudo-gap behavior is observed at relatively high 
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FIG. 1: Transition temperatures T c (thick solid lines) and Td 
(thin solid lines), and zero temperature order parameters 4>o 
(thick dotted lines) and -Do (thin dotted lines) as functions 
of doping level i in a mono-layer system for (a) q = 0.1, (b) 
q = 0.6, (c) q = 1 and (d) q — 1.5. Note that T c and ipo 
coincides wherever To = according to the convention used 
in the text, and they coincides in (b) at all doping levels. 

temperatures, but the associated energy scale extrapo- 
lates to zero near x — 0.19.@ We regard this as yet 
another successful aspect of this mean field theory. 

We now consider a double-layer system, where the only 
modification to the theory is the inter-layer tunnelling, 

/= Y^([ a {x)+f3T^ n \ 2 + ^ n \ i 

n=l ^ 

+[a(x) + qbT 2 ]Dl + h -Di + gD^A 

-J(V^ 2 + c.c), (6) 

where ip n and D n represents the order parameters on the 
ro-th layer, and J = 0.3 is the tunnelling energy scale. 
By symmetry the equilibrium values of the order param- 
eters are independent of the layer index. Therefore the 
analysis can be proceeded in very much the same way 
as for the mono-layer case. The effect of the inter-layer 
tunnelling can be included into a modification of the up- 
per doping limit of the superconducting order, so that 
a(x) — J — > 10(x — 0.33). The discussion then follows 
closely the case of mono-layer. Figs. 2 present the zero 
temperature order parameters and the transition temper- 
atures for the corresponding values of q used in Figs.l. 
Because of interlayer tunnelling that enhances supercon- 
ducting order, the onset doping for the zero temperature 
DDW order is reduced slightly to x = 0.188. Except from 
such details, we find that the qualitative behavior is the 
same as in a mono-layer. In particular bq — g = remains 
to be the special case. 



FIG. 2: The same plot as Figs.l, but for a double-layer system 
with inter-layer tunnelling, because of which the upper doping 
limit of the superconducting order shifts to x — 0.33, and the 
onset of zero-temperature DDW order occurs at x — 0.188. 

Given the success of the theory in dealing with mono- 
and double-layer systems, generalization to more layers 
seems natural and reasonable. The free energy density is 
now written as, 

N ( 8 
f= E W^)+/3T 2 ]|^| 2 + ||V„| 4 

n=l ^ 

-J(CVwi +cc) 

+ [a(x n ) +qbT 2 ]D 2 n + b -D±+ g£> 2 hM 2 ) , (7) 

where N is the total number of layers, and x n is the dop- 
ing level on the n-th layer, which is different from the 
mean doping level x because of the charge redistribution 
effect. || The hole distribution can be roughly described 
by x n = [1 — e/(N — 2)]x for the inner layers 1 < n < N 
and x n = [1 + e/2]x for the outer layers (n — 1 and 
n = iV).@, @ Here e = 0.085, 0.39 and 0.61 for the 
N = 3,4, 5 respectively. We note that the hole distribu- 
tion was measured experimentally mainly in the doping 
regime x > 0.18.|8| Anticipating that e does not change 
appreciably with doping (for a fixed N), we shall extrap- 
olate to obtain the hole distribution in the under-doped 
regime. With the present finite temperature theory, we 
can calculate T c directly from the secular problem in the 
linearized Landau equations, 

[a(x n ) + PT 2 + gD 2 n ]^j n - J(^ n+1 + = 0,(8) 

where D 2 = max{0, — [a(x n ) + qbT 2 ]/b} and n = 
1, 2, N. The condition for the existence of a nontrivial 
solution to the linear homogeneous Eqs.® gives uniquely 
T = T c . We seek symmetric solutions for tp n as a func- 
tion of n, as this yields the highest T c . This reduces 
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the number of independent variables tp„, so that i 
N = 5 we only have to deal with a 3 x 3 matrix d 
minant. The algebra is straightforward. The nume 
result of T c as a function of N is presented in Figs, 
the cases of (a) x = 0.14 (under-doped), (b) x = 
(optimally doped) and (c) x = 0.25 (over-doped) 
each case we used q = 0.6 — 1, as argued to be re; 
able previously. In Fig. 3 (a) T c peaks broadly at iV 
for q = 1. We have checked that for even lower 
ing levels T C (N) increases monotonically (and event 
saturates), to which we shall return shortly. Homolc 
serious in this doping region is not available to the 
of our knowledge but is desirable to check the predic 
In Fig.3(b) T c peaks at N = 3 (4) for q = 0.6J 
and explains nicely the data reported in Ref.0, @ 
Fig. 3(c) T c increases monotonically again. This is 
consistent with experimental data. Indeed, T c (3) = 
and T c (4) = 117iT at x « 0.25 in Table 1 of ref.0. 
perconductors with N > 5 and x > 0.25 are not avai 
in Ref.jl] . 

The non-monotonic dependence in T C (N) is under- 
stood as a cooperative effect from inter-layer tunnelling, 
charge redistribution, and the competing order. For a 
uniform hole distribution, the effect of inter-layer tun- 
nelling is to enhance superconductivity so that T c in- 
creases with the average layer-coordination number 2 — 
2/N (as has been checked but not shown here). This ac- 
counts for the initial enhancement in T C (N) for all doping 
levels. In the very under-doped region the hole density 
is low despite of the charge redistribution effect. In such 
cases, the competing DDW order is robust and does not 
change appreciably. Therefore, apart from a global sup- 
pression by DDW, T c should increase with N. This is 
the case in Fig. 3(a) if one takes q = 0.6 as the appropri- 
ate parameter. In the very over-doped region, the hole 
density is so high that DDW order is absent (at least 
up to T c ), and only the inter-layer tunnelling matters, 
yielding an increasing T C (N) again, as in Fig. 3(c). The 
situation is quite different near the optimal doping. At 
the average doping, DDW is weak or absent. But with 
increasing N, charge redistribution introduces hole-poor 
inner layers with DDW order, counteracting the effect of 
inter-layer tunnelling and eventually suppressing T c , as 
in Fig.3(b). 

To summarize, we proposed a finite temperature Lan- 
dau theory that describes competing orders and inter- 
layer tunneling in cuprate superconductors as an im- 
portant extension to a corresponding theory at zero 
temperature. For given parameters, our theory deter- 
mines T c without any ambiguity. In mono- and double- 
layer systems we discuss the relation between zero tem- 
perature order parameter and the associated transition 
temperature in the presence of competing orders, and 
discuss the puzzling experimental fact that the pseudo- 
gap temperature is much higher than the corresponding 
energy scale near optimum doping. Applying the the- 
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FIG. 3: The layer-number N dependence of the superconduct- 
ing transition temperature T c . (a) x — 0.14, (b) x = 0.20, (c) 
x = 0.25. 

ory to multi-layer systems, we calculate the layer-number 
dependence of T c . In a reasonable parameter space the 
result turns out to be in agreement with experiments. 

We should emphasize that just as in Ref. 2] DDW 
is used as an example of the competing order 
to superconductivity. The discussion in this pa- 
per holds for any other forms of competing orders. 
However, DDW has the special property that it 
yields a band-gap with d-wave symmetry in the mo- 
mentum space, 4] in accordance with angle-resolved- 
photoemission measurement. 9] While this seems to be 
appealing, other consequences that follow the DDW pic- 
ture are still in debates. For example, the temperature 
derivative of the superfluid density at low temperatures 
was claimed to be insensitive to the doping deep in the 
under-doped limit of YBCO thin films. |l0j but another 
group reported that in YBCO crystals it tends to di- 
verge with under-doping. ^ While DDW picture is con- 
sistent with the latter, it can not account for the former 
results. [l2j The fate of the DDW picture depends on fur- 
ther experimental clarification. 
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